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Abstract: The cometary tray of the NASA Stardust spacecraft’s aerogel collector has 
been examined to study the dust that was captured during the 2004 fly by of comet 
81PA¥ild-2. An optical scan of the entire collector surface revealed 256 impact 
features in the aerogel (width > 100 pm). 20 aerogel blocks (out of a total of 132) 
were removed from the collector tray for a higher resolution optical scan and 186 
tracks were observed (track length > 50 pm and width > 8 pm). The impact features 
were classified into three types based on their morphology. Laboratory calibrations 
were conducted which reproduce all three types. This work suggests that the cometary 
dust consisted of some cohesive, relatively strong particles as well as particles with a 
more friable or low cohesion matrix containing smaller strong grains. The calibrations 
also permitted a particle size distribution to be estimated for the cometary dust. We 
estimate that approximately 1200 particles bigger than 1 pm struck the aerogel. The 
cumulative size distribution of the captured particles was obtained and compared with 
observations made by active dust detectors during the encounter. At large sizes (>20 
pm) all measures of the dust are compatible, but at micrometer scales and smaller 
discrepancies exist between the various measurement systems which may reflect 
structure in the dust flux (streams, clusters etc.) along with some possible instrument 
effects. 
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1. Introduction 

The use of silica aerogel (originally developed as a low density material by Kistler 
1931, with ultra low densities available in recent decades), as an ultra low density 
transparent solid to capture relatively intact samples of projectiles impacting at speeds 
in excess of a few km s' 1 was demonstrated in the laboratory (Tsou et al. 1988) and a 
recent review of subsequent developments is given by Burchell et al. (2006). Shortly 
summarised, silica aerogel is a solid foam, in effect a dried SiOi gel with an open 
network of pores. The low densities obtainable today for aerogel (which range from 
just over 1 to approximately 500 kg m 3 ), are the consequence of this pore space 
inside an otherwise solid material. Projectiles travelling at hypervelocity speed (of 
order km s’ 1 ) normally undergo severe damage on impact with a solid (non-porous) 
target of density much greater than aerogel. Typically, during such impacts the front 
of the projectile slows but is then crushed by the rear of the projectile (which has not 
yet slowed) and the target material cannot initially move away from the impact point 
to respond to the impact, thus generating a shock wave. The result is that extremely 
high densities and shock pressures are generated in both the projectile and target 
materials. During adiabatic release from this shocked state (which can involve 
pressures of tens or even hundreds of GPa), heating of the materials occurs. Most of 
the projectile is then vaporised, and only a few % of it can be subsequently found at 
the impact site, usually in the form of a melt lining part of the crater which forms in 
the target. 

In aerogel however, the impact process is different and is more of a 
“penetration” event than a cratering one. The peak pressures generated upon collision 
with such a low density medium are insufficient to melt typical silicate-impactors, 
sometimes even insufficient to cause projectile fragmentation. Such projectiles tunnel 
into the aerogel, losing speed as a continual process (see Anderson and Ahrens 1994 
for a general model for capture of hypervelocity particles in mesoporous foams and 
Dominguez et al. 2004 for a recent model of capture specific to aerogel). At some 6 
km s’ 1 , the projectile may penetrate a factor of 100 - 200 deeper into aerogel 
compared to crater depth in solids such as rock or metal, producing a long, readily 
visible penetration track with a captured and typically un-melted particle at the end. 
During capture, there is some heating of the aerogel which can melt and form a 
molten wrap around the particle and some of the projectile surface is heated and shed 
as an ablative process. Depending on the nature of the impactor, afterwards a 
relatively intact particle can be found at the end of a track in the aerogel. How much 
of the particle is captured depends on impact speed, aerogel density and projectile 
structure and composition, but at 5 - 6 km s’ 1 as much as 60, 70 or even 100% of 
relatively competent particles some tens of pm across can be captured (see Burchell et 
al. 2006). Alternatively, the impact process can so severely disrupt the particle that it 
effectively breaks apart into a large number of fragments lining the wall of a wide 
cavity in the aerogel. And some inhomogeneous particles give a combination of these 
two capture possibilities, resulting in tracks which are a mixture of the large, broad 
cavity plus discrete subsidiary tracks emerging from it. This is discussed in more 
detail later. 

Given that silica aerogel is highly transparent, and that the particle fragments 
(large or small) are easily observable and readily extractable for detailed 
characterization with state-of-the-art analytical instruments, aerogel became the 
collector medium of choice for retrieving cosmic dust, i.e. dust in space which is 
travelling at speeds of km s’ 1 relative to the observer. Several space missions have 
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deployed aerogel dust collectors in Low Earth Orbit (again reviewed in Burchell et al. 
2006). 

The NASA Stardust mission to comet 8 1 P/Wild-2 deployed an aerogel dust 
collector (Brownlee et al. 2003) as it flew past the comet at 1.86 AU from the Sun, 
with a relative encounter speed of 6.1 km s' 1 (Brownlee et al. 2004; Tsou et al. 2004). 
At this heliocentric distance the comet was actively ejecting dust and gas. The speed 
of the dust grains (with respect to the comet) was low compared to the fly by speed of 
the spacecraft so the impact on the aerogel was at a constant encounter speed of 6.1 
km s" 1 and normal incidence. An active impact sensor on the spacecraft’s leading edge 
(the Dust Flux Monitor Instrument - DFMI) indicated that during the encounter dust 
grains were indeed striking the spacecraft and it predicted a flux for the entire aerogel 
collector of 2800+500 cometary grains larger than 15 pm in diameter (Tuzzolino et al. 
2004), although an impact ionisation detector (sensitive to smaller particles) indicated 
a much lower flux (Kissel et al. 2004) 

The successful return to Earth of the Stardust aerogel collector in January 
2006 provided access to the cometary material captured in the aerogel. The first 
reports (Brownlee et al. 2006; Horz et al. 2006 and other papers in the same issue of 
Science) have already provided a wealth of detail from the Preliminary Examination 
(PE) post flight analysis effort. In this paper we provide an in-depth analysis of the 
characteristics of the tracks observed in the Stardust cometary aerogel, combined with 
a description of experimental impacts that simulated the Stardust environment, thus 
providing insights and “calibrations” for the Stardust observations. 

2. Laboratory Calibration 

2.1 Laboratory Impacts With Glass Beads 

The calibration of track size as a function of impactor size or mass is greatly 
aided by the fact that the encounter speed of Stardust was a constant 6.1 km s’ 1 at 
normal incidence to the collector surface. The experiments employed the same soda- 
lime glass beads that were used to calibrate the crater size in Stardust aluminium foils 
(Kearsley et al. 2006). These glass beads were supplied by Whitehouse Scientific 
Ltd., and their characteristics are given in Table 1. In Kearsley et al. (2006) particle 
sizes were measured independently of the supplier’s nominal values and good 
agreement was found. Using the University of Kent two stage light gas gun (Burchell 
et al. 1999a), samples beads in three size ranges were fired into Stardust grade aerogel 
manufactured as part of the same batches used in the Stardust cometary aerogel 
collector. This aerogel did not have a uniform density, instead it varied across the 3 
cm thickness of the blocks, from 5 kg nT at the front face to 50 kg m" at the rear face 
(Tsou et al. 2003). The impact speeds were measured for each shot to better than 1% 
and given in Table 1. Tracks were observed with the classic carrot shape, i.e. entrance 
hole widening quickly to a maximum track width just below the surface and then 
tapering with (non-parallel) straight-line sides to a near point where the particle was 
captured (see schematic in Figure 1). Note that the maximum track width is not at the 
entrance hole, but occurs at some small depth below it. 

For each size of projectile several tracks were measured under the microscope. 
The various parameters measured are defined as follows: Track length is from the 
straight line distance from centre of the entrance hole in the surface plane of the 
aerogel to the centre of the captured particle. The entrance hole diameter is the 
average of several measurements; the holes are not necessarily circular so two 
orthogonal diameters are taken and averaged, any narrow fractures in the surface 
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plane which emanate from the main hole are ignored. The maximum track width is 
measured transverse to the main track axis (i.e. track length) at its widest point and 
the depth of this position below the original aerogel surface is also recorded. The 
captured particle may not be spherical; accordingly two orthogonal diameters are 
taken and averaged. The average values for the measurements of each quantity are 
given in Table 2. The uncertainty given on each average value is the la value of the 
distribution in each case. On average the maximum track width occurred 13% of the 
total track length below the entrance hole (although there is a large scatter on this). 
The data for track length (L), entrance hole diameter (EHD) and greatest width (W) 
are plotted vs. original particle diameter (OPD) in Fig. 2 (various fits are shown and 
are described as follows). In Figure 2a, it is found that 

F = -(1287+669) + (269+19) x OPD, (1) 

where units are in pm. This can be re-arranged to obtain: 

OPD = [L +(1287+669)] / (269+19), (2) 

where all units are pm. If the fit is forced to pass through the origin, then we obtain 
that F = (269+12) x OPD (dashed line in the figure), very similar to the original fit. 

The data for EHD’s are shown in Figure 2b, along with three fits. The first fit 
is a linear fit which gives 

" EHD = (14+46) + (5.0+1. 4) x OPD. (3) 

Again units are pm and the results can be rearranged as: 

OPD = [EHD- (14+46)] /(5.0+1.4). (4) 

If the linear fit is constrained to pass through the origin we obtain EHD = (5. 3+0. 6) x 
OPD, shown as a dashed line in the figure and very similar to the original linear fit. 
The third fit considered in Fig. 2b is a polynomial fit. As can be seen in the figure, 
although in the size range used it tracks the data points themselves, this fit does not 
approach the origin at small values of OPD, a non-physical result for small sized 
particles. This is typical of fits to any of the calibration data, where a linear fit is 
always found to well describe the data. With this caveat about failing for small 
particles, the fit yields 

EHD = (60+71) + 1.38 x OPD + 0.057 x OPD 2 , (5) 

(units are pm). This can be solved as a quadratic to give the original particle diameter 
for a particular value of entrance hole diameter (EHD). 

OPD = {-1.38 + V[ 1 .904 - 0.228x(60 - EHD)]} / 0.114 (6) 

In Fig. 2c the greatest track width data are shown and are again fitted twice. 
Both fits are linear. The first (solid line) only uses the 3 data points and gives 

W = -(36+50) + (13.3+2.3) x OPD, (7) 

units are in pm. Rearranging gives 

OPD^= [(36+50) + W] / (13.3+2.3). (8) 

The intercept on this fit is a finite value but is within la of zero. Accordingly, we refit 
(dashed line) requiring the linear fit to pass through the origin and obtain 

V = (12.1+0.5) x OPD, (9) 

which rearranged yields 

OPD = W/ (12.1+0.5), (10) 

where units are pm. 

As well as the directly measured quantities, it was also possible to estimate the 
volume of each track. This was done by approximating the track shape as a series of 
congruent frustra. A frustum is a symmetrical truncated cone as shown in Figure 3a. 
Two such shapes can be combined to approximate a carrot track shape (Fig 3b). The 
key measurements of a frustum were all available from the track measurements. 
These were then used to obtain the volume (V) of the 3 -dimensional track. On 
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selected tracks volume measurements were made by slicing the track into 20 or 40 
individual cells and finding the volume of each. Agreement in the results with the 2 
frustrum method was at the level of 10 - 15%, indicating reasonable robustness in the 
method. The average results for volume are given in Table 2 and plotted vs. impact 
kinetic energy (KE) in Figure 4. The KE’s were obtained for each sample using the 
nominal particle size and the mean mass, combined with the speed of each shot. This 
yields energies of 3.5 x 10’ 5 J, 9.1 x 10" 4 J and 6.02 x 10" 3 J respectively for the 11.58, 
35 and 63.8 pm diameter samples. The fit to Figure 4 gives 

V = -(0.011 ±0.009) + (600+170) xKE, (11) 

where units are mm for volume and J for KE. 

Rearranging this gives 

KE = {V + (0.011 +0.009)} / (600+170). (12) 

where again units are mm' for volume and J for KE. This suggests 1.7 mJ of energy 
are needed to excavate each cubic mm of aerogel, with an accuracy of about 28%. 
That the intercept is only just over la away from zero, suggests this relation also 
holds for smaller tracks than those measured here. 

In general there is a degree of surprise about these results. For example it has 
been shown previously (e.g., Horz et al. 1998; Burchell et al. 2006) that for impacts at 
equal speed by similar sized particles, track length increases as aerogel density 
decreases. Here, impacts of different sized projectiles into density gradient aerogel 
have produced a simple dependence of, for example, track length on particle size. 
Detailed modelling is underway to evaluate the origin of these dependences and to 
determine the actual density gradient in the top 10 mm of Stardust aerogel tiles. It 
should also be noted that with only 3 data points in each calibration data set here, 
although linear fits are shown to be good descriptions of the data, any fine non-linear 
dependencies are not discernible rather than completely excluded. 

Using these results it is possible to measure track length, entrance hole 
diameter, or greatest track width and obtain an estimate of the diameter of the original 
impacting particle on the Stardust collector. In the case of track volume the calibration 
gives an impact energy (based on the volume of aerogel excavated). For the known 
impact speed of Stardust, a particle mass can then be found. If a density is assumed 
the particle diameter can then be obtained. 

The above relationships were obtained for soda-lime glass beads (density 2.4 g 

a 

cm" ) impacting the aerogel. Whether these can then be directly applied to other types 
of particles requires consideration. In the case of entrance hole diameters, it is 
reasonable to assume that the critical parameter is particle diameter. For example, 
Burchell et al. (1999b) have shown that 100 micron diameter particles of olivine, soda 
glass and iron all leave similarly sized entrance holes in aerogel. This seems 
reasonable if we consider the analogy as penetration of a very thin sheet; at 
hypervelocities the hole punched in thin films gives a good measure of particle cross- 
sectional area (Horz et al. 1994). If taken further, this analogy may suggest that 
similar to the thin film, the continuous aerogel medium does not exhibit a major flow 
field around the surface entrance hole. Track length should not however be taken as a 
good indictor of particle size, as this may well depend on additional properties of the 
particle (e.g. density, which will affect impact energy or melting point that will affect 
the mass loss during penetration). Greatest track width is also not an ideal parameter, 
although it is accurate for tracks similar in shape to those obtained in the calibration, 
in the Stardust aerogel other track morphologies are also seen. Track volume may be a 
good indicator of impact kinetic energy as the data suggest a constant excavation 

a 

energy per mm of aerogel for the soda lime glass impactors used in the present 
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experiments. This assumes no extra source of energy (i.e., from chemical 
decomposition of the projectile during capture). This will be returned to later when 
applied to real data from Stardust. 

2.2 Laboratory Impacts With a Variety of Minerals and Basalt. 

As well as glass beads, impacts on Stardust cometary grade aerogels were obtained in 
the laboratory with 4 types of natural materials. These were to simulate impacts from 
a typical range of materials possibly representative of those likely to be found at a 
comet. Details of the materials are given in Table 3. The olivine and basalt grains 
gave carrot tracks in the aerogel similar to Fig. 1, these are called Type A. The 
pyrrhotite grains gave Type B tracks, i.e. tracks with a slightly broader initial cavity 
which then tapered down to a carrot track. This has the appearance of a champagne 
flute wine glass without the base (i.e. only a bowl and a stem region). This is shown in 
Fig. 5a. The lizardite grains however, gave Type C tracks with a broad cavity with no 
individual tracks emerging from it (Fig. 5b); the walls of these cavities were lined 
with fine fragments of particles, some of which had penetrated slightly into the 
aerogel beyond the edge of the cavity. Several such broad cavity tracks were seen in 
the aerogel from the lizardite impacts. In Fig. 5b a typical cavity is shown at cm size 
scale, we also observed similar impact features at mm scales. Thus this behaviour is 
not size dependent. However, in the lizardite shot there were also some impacts which 
gave Type B tracks with discrete fragments (of lizardite) at the end of the stylus. The 
size distribution of these trapped grains (and the associated large scatter in track 
length) did not reflect the size of the initial particles, suggesting they may have 
broken up during launch before impacting the aerogel. 

There are thus three types of tracks observed in these shots and the 
morphological classification used is that which was also used in Horz et al. 2006. 
Type A (soda lime glass, basalt, olivine) is the classic “carrot” shaped aerogel track, 
which initially broadens to a maximum width below the initial impact point, then has 
(non-parallel) near straight lined walls which come to a point (with the trapped 
particle nearby). Type B (pyrrhotite and some lizardite) are also relatively slender 
tracks, but have a broader initial cavity than the glass projectiles and curved walls 
which taper to a thinner stem from which emerges a carrot- shaped track, terminating 
near a captured particle. Type C tracks (some lizardite particles) have only a broad 
cavity and no stem emerging from it, with no large discrete fragments of the initial 
particle beneath the cavity, instead multiple small fragments line the cavity. 

There has been previous experimental work which has produced aerogel tracks 
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of various shapes. Horz et al. (1998) observed impacts in aerogel (density 20 kg m' ) 
in the laboratory at speeds up to 6 km s' 1 . As well as carrot shaped tracks from glass 
beads, they also observed Type B and C tracks from impacts of compressed dry cocoa 
powder (made from micron sized grains). Type B tracks were obtained by mixing 
aluminium or glass spheres (50 pm) with the cocoa power, and Type C by pure cocoa 
powder projectiles. When examined afterwards, the cocoa powder was found to line 
the wall of the cavities. Their conclusion was that low cohesion projectiles would 
leave bulbous, Type C cavities, but cohesive projectiles produce carrot-shaped Type 
A tracks: mixtures of the two types of material inside one inhomogeneous particle 
would produce a track type which was a cross between these two extremes, i.e. a Type 
B track. There is thus good agreement with the types of tracks found in the mineral 
shots here in Stardust grade aerogels and the work of Horz et al. One difference is that 
the bulbous cavities in the lizardite shots here may arise from a high volatile content 
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of the projectile (lizardite is typically 14% water by mass) whereas in Horz et al. it 
was due to low cohesive strength of a fine grained matrix. At potentially higher, albeit 
unknown, impact speeds, in low Earth orbit Horz et al. (2000) also observed similar 
track morphologies in aerogel exposed on the Mir Space Station. Separately, 
Kitazawa et al. (1999) proposed a scheme for track categorisation. For impacts at > 4 
km s' 1 , they found three track types (similar to those here) from impacts into aerogel 
of density 30 kg m 3 . There is thus good agreement between three separate sets of 
experiments indicating that three major types of tracks are obtained from impacts in 
aerogel. In the remainder of this work we use the names Type A, B or C to label each 
type. The transition between types (and variation inside types) is totally gradational. 

Although three types of tracks are found, the entrance holes appear similar (a 
typical entrance hole is shown in Figure 6). Thus entrance hole does not differentiate 
between track type. Although the mineral grains were not monodispersive, the range 
of entrance hole sizes in each shot is compatible with that expected from the known 
particle size range and the calibration for soda lime glass (Fig. 2b and eqn. 3 and 5). 
This indicates that entrance hole diameter is probably the most suitable parameter to 
use for determining particle size from impacts in aerogel of unknown particles. 
However, due to the highly transparent nature of aerogel and some intrinsic surface 
roughness, it can be difficult to focus on the surface face of an aerogel block and 
obtain a well resolved image of an entrance hole, especially at small track sizes. 


3. Stardust Cometary Tray Tracks 
3.1 Level 2 scan 

The returned Stardust aerogel samples were examined in a dedicated Class 100 clean 
room at the NASA Curatorial Facilities in Houston. The cometary aerogel tray 
contained 132 cells of aerogel, 130 of which were rectangular with slightly rounded 
comers (surface of nearly 4 cm x 2 cm) but two were trapezoid shaped, with width 2 
cm and on the long sizes lengths of 3 and 2 cm. The total surface area of aerogel was 
1039 cm". All cells were 3 cm deep. The aerogel density was not less than 5 kg m' at 
the front surface, increasing to no more than 50 kg m 3 at the rear. Full details are in 
Tsou et al. (2003). The aerogel blocks were held in an aluminium frame with an 
appearance like a tennis racket (see Fig 7), with each cell held in place by soft 
aluminium foil along the walls of the cell and over the sides of the supporting frame. 
The blocks could be removed by cutting these foils on the flat top of the frame and 
pulling the foil from the rear to slide the block out. The exposed surfaces of these foils 
were also analysed for impact features (Horz et al. 2006; Kearsley et al. 2007). 

Prior to any sample processing, all cometary aerogel surfaces were imaged, 
via optical microscope/CCD at xl6 (and in some cases also at x20) magnification, 
typically producing some 42 images for each individual aerogel block at a resolution 
of approximately 3.5 (or 2.7 respectively) pm/pixel. This is referred to as Level 2 
Stardust Photography documentation. Because the aerogel cells still resided in the 
modular openings of the collector tray, this photography depicts all tracks in plan 
view (from above or at 15° inclination). Systematic analysis of these mosaics and of 
individual framelets was conducted for all 132 aerogel blocks, aided by the iterative 
microscopic inspection of actual tiles to clarify the nature of some features that 
seemed ambiguous in the images. An example entrance hole is shown in Fig 8 (note 
that tracks are numbered in the style Cnnn-Tm, where nnn is the 3 digit number 
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referring to the aerogel piece in the collector tray, see Fig. 7, and m is an integer 1, 2, 
3.... for the first, second, third, etc., track found in that block). Although the imagery 
was focussed on a cell’s surface, determination of the entry hole was difficult in most 
cases, as features at some depth in the transparent medium often interfered. Thus the 
parameter most readily measured in a systematic fashion in plan view is the maximum 
diameter (W) of individual tracks, modestly below the surface, yet still in reasonable 
focus. The threshold for reliable track recognition during this tray wide survey was 
empirically set at W = 100 pm, a compromise between actual resolution, surface 
roughness, and operator time. This survey produced a total of 256 tracks > 100 pm 
diameter; Fig. 9 illustrates their size frequency. The accuracy of the measurements 
was ±10 pm, and accordingly at the smallest sizes the data are grouped in 10 micron 
bins. The greatest track width recorded was 9.961 mm. The distribution of these 
features across the surface of the aerogel collector is shown in Figure 10. Given 132 
aerogel blocks and 256 features, approximately 2 per block are expected for a random 
distribution. 

To generate a preliminary impact flux for the whole cometary aerogel 
collector a calibration has to be applied. The relation of eqn. 8 is used to provide 
estimates of projectile diameter Dp (note that for greatest widths of diameter >100 
pm there is little difference between use of eqn. 8 or 10). The resulting cumulative 
particle size was given in Horz et al. (2006) and is shown in Fig. 1 1 . As expected the 
distribution cuts off at about 10 pm in D P and extends up to almost 1 mm. Ignoring 
the largest few particles (where small statistics may cause fluctuations), a fit to the 
data from 10 to 200 pm gives 

Cumulative size distribution per m 2 = 47100 Dp' 122±0 06 eqn (13) 

At larger sizes the slope in the cumulative flux increases somewhat, but as 
stated this may be due to the statistics of small numbers of events. The lack of a roll 
off in the flux at smaller sizes indicates a fairly complete degree of scanning. Based 
on this flux we obtain that 180 particles greater than 15 pm in diameter struck the 
aerogel. 

However, some caution is required in fully accepting this estimate of the flux. 
The calibration used was obtained from the experimental Type A tracks. However, as 
will be seen below, the data here are a mixture of track Types. As already stated it is 
not certain if maximum track width for a given value of D P is independent of track 
Type. Type B tracks are wider than Type A, by perhaps a factor of 2 (based on rough 
estimates made from the calibration shots of lizardite). So a fraction of tracks in the 
cumulative flux in Fig. 1 1 may have had their D P overestimated by a similar factor. 
Since Type A tracks dominate at small track lengths the true slope in Fig. 11 is 
probably slightly steeper than that in eqn. 13. 

3.2 Level 3 Scan 

After the whole tray survey, 20 cometary cells of aerogel were initially removed from 
the tray and their analysis is given here (shown as shaded blocks in Fig. 7 and listed in 
Table 4). (Note: additional cells have been removed subsequently either for analysis 
or long term curatorial storage). The removed cells were then imaged side on, where 
each tile was rotated some 15° to permit viewing of the exposed surface from below, 
which substantially improved -compared to perfectly orthogonal views- the 
recognition of very small tracks. This so called Level 3 Stardust Photography 
documentation was conducted at magnifications of up to xl6, x25, x32, x50, x63 or 
even xlOO, with a resolution of 3.4, 2.2, 1.7, 1.1, 0.88 or 0.55 pm per pixel 
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respectively, modestly higher than Level 2. Also, each tile was photographed in 
discrete depth intervals, progressing in 2.5 mm steps from the front to the rear side, 
resulting in 8 “slices”, each slice consisting of some 64 images. The purpose of this 
photography was to produce a permanent record of the size and X/Y location of the 
track population in a given tile. Additionally, most of the large tracks were 
photographed individually, with optimally adjusted focal plane. The extracted blocks 
represented some 15% of the total aerogel in the cometary tray, sufficient to provide a 
high resolution survey of tracks in the collector. The rest of the tray is currently being 
preserved by NASA as a resource for future analysis. 

A total of 206 possible track-like features were initially noted in the Level 3 
images. One feature was excluded as it had impacted on the edge of the supporting 
aluminium frame and the resulting debris had spread into the aerogel and 19 features 
were discarded because they were somewhat ambiguous and not verified by a second 
observer. This gave a set of 186 tracks in the 20 aerogel blocks. Due to their small 
size and faint images, 9 of these tracks could only have their dimensions measured to 
worse than 20% accuracy. The remaining 177 tracks were measured to higher 
accuracy, with the largest measured to better than 1%. The minimum track length was 
55 pm and the smallest maximum track width measured was 6 pm. 

The tracks were visually categorised into Type using the 186 track data set. 
Examples are shown in Fig 12 (Type A), Fig. 13 (Type B) and Fig 14 (Type C). In 
Figure 12b the captured particle at the end of a Type A track is shown. In Figure 13, 
very thin tracks can be seen emerging almost radially (with a slight downward tilt) 
from the main cavity. Fine grains of particulate material line the main cavity, and are 
seen along these fibre like tracks and the stylus, as well as there being a terminal 
particle at the end of the stylus. Many of these fragments are micron scale and suitable 
for analysis by a variety of techniques. Thus a single Type B track may contain many 
identifiable particle fragments suitable for compositional studies. There is a very fine 
division between Types B and C, several tracks listed as B are close to type C, 
possessing only very short styli emerging beneath the main cavity. 

Averaged over all sizes, 65% of tracks were Type A, 33% Type B and 2% 
Type C. However, this was found to be size dependent, and a breakdown of the 
percentage of Type A and B vs. track length is given in Fig 15 and given numerically 
in Table 5. For small track lengths (< 100 pm) all tracks are given as type A. There is 
a problem here as the resolution is starting to become comparable to track width, and 
although the tracks do not appear to have a transition region along their length (which 
would indicate a Type B track if it were seen), higher resolution imaging is required 
to determine if this is indeed the case. For tracks greater than 500 pm length, 50% are 
A’s and 50% B’s. In the range 100 to 500 pm track length, there is a gradual 
transition from one regime to the other. At larger track lengths (>1500 pm) the small 
statistics prevent meaningful comparison. The 4 Type C tracks were spread over all 
lengths in proportion to where the bulk of the tracks occurred, indicating no favoured 
size scale. 

In some impacts of Type A the incident particle appears to have split during 
capture, and resulted in more than 1 terminal track (or “stylus”) per impact. This is 
also true for Type B tracks where several styli can be seen emerging beneath the main 
cavity (and are clearly distinct in size from the radial hair like features seen in Figure 
13). The distribution of how many styli are found in each track Type is shown in Fig. 
16. There is a clear difference between track Types. Type A tracks are dominated by 
single terminal particles (86%). Type B tracks have one terminal particle just under 
50% of the time, two in 1/3 of occurrences and 21% of Type B tracks have 3 or 
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more terminal particles. This suggests a different nature to the particles causing Type 
A and B tracks. Type A particles are stronger materials, well consolidated into usually 
a single grain. Type B particles are weaker in some fashion, either more loosely 
bound or more volatile rich, which are more readily disrupted into fragments during 
capture. 

Each track was then individually measured. The key measurements were: 
Total track length (L) from entrance hole to the deepest penetrating discrete fragment 
(or bottom of cavity for Type C tracks) and entrance hole diameter (EHD) and 
maximum track width (W). In addition, for Type B and C tracks, the main cavity was 
divided into several contiguous frustra and the necessary heights and widths obtained. 
For all tracks, if more than 1 stylus was observed each was measured separately (and 
treated as being similar to a “carrot” shaped track, i.e. it had a base diameter and 
length which served to define a cone shape for the stylus). From the set of 
measurements made for each track, an estimate of track volume (V) was then 
obtained, treating each track as a set of frustums or cones as required (similar to the 
treatment of the calibration data). 

The distribution of maximum track widths vs. lengths is shown in Fig. 17. The 
track lengths ranged from 50 to 10,000 pm. Track widths ranged from 4 to 2000 pm. 
A strong correlation was observed between track width and length (indicated by the 
near diagonal trend in the data on Fig. 17). The data displayed in Fig. 17 have been 
grouped by track Type. A division between the types is possible (with a few 
exceptions close to the boundaries). The division between Types A and B is given by 

W = 0.065 xF 11 , eqn. (14) 

and between Types B and C by 

W = 0.20 x F 1 ' 1 . eqn. (15) 

As note above, at short track lengths (F < 100 pm) only Type A tracks are observed. 
It is possible that for the shortest tracks, the imaging resolution is insufficient to 
clearly distinguish between Types A and B. The tracks appear at current resolution as 
“fat” cylinders with only some slight tapering towards their terminus. They are thus 
not apparently Type B, but are however relatively broader than the Type A track seen 
at larger scales. This assignment should thus be taken with caution and will be the 
subject of future work at higher resolutions. 

The ratio W/F is informative for tracks, as it indicates how relatively thin or 
fat a track is. W/F is shown vs. track length in Fig. 18. The lines shown on Fig. 18 
divide the date into distinct regions. At low track lengths (< 100 pm, shown by a 
vertical solid line) only type A tracks are seen (as noted earlier). The main population 
of Type A and B tracks are mostly separated by a horizontal line shown at W/F = 
0.11, whilst the main population of Type B and C tracks are separated by a horizontal 
line shown at W/F = 0.35. Ignoring the region F < 100 pm, track Type A thus has 
W/F <0.11, Type B has 0.11 < W/F < 0.35 and Type C has W/F > 0.35 . Some Type 
B tracks do intrude into the region defined for Type C. Renewed examination of these 
shows that the terminal track in these Type B tracks is visible, but is typically very 
short and narrow (in some images similar short tracks are seen radiating near laterally 
from higher up the main cavity), indicating the presence of a very small fragment of 
the original particle. 

In Fig. 19 the track volumes are shown vs. track length for all tracks. In Fig 
19a are Type A tracks, and the data were fitted to yield 

V = 0.043 x l ( 2 47±0 ° 8) , eqn. (16) 
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where L is in pm and V in pm . The behaviour of the data for L < 100 (am seems to 
differ somewhat for that at larger track lengths. Accordingly, this data was excluded 
from the fit and a new relation obtained for Type A tracks (L > 100 pm) of 

V = 0.012 x l < 2 ' 66±0 09) . eqn. (17) 

As might be expected, the power in eqn. 17 has increased slightly compared to eqn. 
16. Both fits are shown on the plot (all data solid line, L > 100 pm as a dashed line). 
In Fig. 19b the data are given for Type B tracks and these yield a fit: 

V = 0.0072 x L (3 01±010) , eqn. (18) 

where again L is in pm and V in pm 3 . In Fig. 19c, all data are shown combined into 
one plot. Ignoring the region for L < 100 pm (boundary marked by a vertical line), we 
find that the boundary between Type A and B tracks is given by 

V = 0.003 xL 3 , eqn. (19) 

and between Types B and C by 

V = 0.05 x L 3 . eqn. (20) 

In most cases attempts to directly measure the size of the captured particles 

from the Level 3 images failed as the particles were at or below the limits of resolving 
and measuring such tiny grains in the images. More detailed images of the captured 
particles were made for a few tracks as part of the compositional analysis studies 
which are reported elsewhere in this volume. 

The measured values for track volume were then combined with the 
calibration data (eqn. 12) to obtain estimates of the pre-impact particle diameters in 
each track. The assumption made was that all particles had a density of 2400 kg m" . 
The results are shown vs. track length in Fig. 20. The projectile sizes ranged from 0.4 
to 100 micron. The frequency of particle size for each track Type is given in Fig. 21. 
The cumulative size plot is shown in Fig. 22. At small sizes (< 1 pm) the cumulative 
size distribution flattens off, indicating a loss in scanning efficiency rather than a lack 
of smaller particles. At large particle diameters, the data are susceptible to the 
statistics of small numbers and the cumulative curve shows a sudden change in slope. 
Finally, the total size distribution of Fig. 24 was fit over the range 1-100 pm, 
yielding 

Cumulative size distribution per nr = 11900 Dp‘ 0 76±0 01 eqn (21) 

The number of particles bigger than 1 pm diameter striking the Stardust 
aerogel tray is predicted from eqn 21 to be approximately 1200 and the number of 
particles larger than 15 micron diameter is predicted to be 158 particles. This latter 
value is compatible with that predicted by the Level 2 scan of the entire tray which 
predicted 180 particles greater than 15 pm diameter. And as shown in Horz et al. 
2006, the size distribution (ignoring clustering which is discussed in the next section 
and can lead to significant variations in flux on small surface areas) is compatible 
with the average size distribution obtained from the craters on the aluminium foils at 
that size scale. This is however, substantially less than the (2800+500) particles 
greater than 15 pm predicted by the DFMI during the fly by in 2004 (Tuzzolino et al. 
2004). 


4. Discussion 

4.1 Calibration and classification of Stardust aerogel tracks 

The calibration is based on two data sets. The first concerns measured track 
parameters (length, volume etc.) obtained from mono-dispersive glass beads fired into 
Stardust grade aerogels in the laboratory. The second uses polydispersive grains of 
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various materials, again fired at Stardust grade aerogels in the laboratory. In both 
cases the data sets are revealing and provide the necessary insights into capture in the 
Stardust aerogels. However, both could be further extended. The present limitation is 
operational and related to the lack of more Stardust grade aerogel. This situation may 
be improved in the near future when the original Stardust aerogel manufacturing 
procedure may be repeated to provide more aerogel for calibration purposes. It will 
then be possible to conduct calibrations with mono-dispersive particles of a range of 
densities as has been done for the cratering in Stardust foils (Kearsley et al., 2007). 
Similarly, fine details of the capture mechanism and how it changes the nature of the 
impactor (e.g. processing of organics during capture) can be studied in more depth. 

The data herein concerning tracks in the Stardust aerogel from the cometary 
collector side were obtained during the PE period up until August 2006. They thus 
represent a first look at the aerogel. Work is on-going, and more data will be available 
later as more blocks are removed for analysis. However, the current data sets are large 
enough to show the main trends in the analysis, although in some cases (e.g. the Type 
C tracks) there are still low statistics which if improved may reveal more. In general it 
is unlikely that many more large tracks will be found. However, with greater statistics 
and particularly higher resolution images for the very short tracks, the track 
morphological classification can almost certainly be improved and the nature of the 
tracks at the smallest scales clarified further. One major point to note however is that 
the data refer only to the cometary collector tray. The interstellar collector tray still 
awaits an equivalent analysis when the first imaging of the tray is complete by the 
Stardust@home analysis effort. 

One notable feature of the observed tracks is their distribution across the 
collector as observed in the Level 3 scan of extracted aerogel cells shown in Fig. 22. 
The observed number counts per cell are not compatible with a random distribution. 
This effect is strongly correlated with feature size (i.e. particle size). For larger 
particles the distribution is nearly compatible with a random distribution, but for the 
smallest particles this is not the case. This is discussed in detail in a separate paper 
(Westphal et al. 2007). 

4.2 Cumulative Size Distribution and Flux 

The cumulative size distribution obtained from the data reported in depth in 
this paper can be compared to that obtained by analysis of the cratering in the 
aluminium foils on Stardust and the flux measured during the cometary encounter by 
electronic instruments read out in real time. All methods can yield a size distribution, 
but the active instruments (DFMI) also provides timing information reflecting on the 
spatial distribution of the dust near the comet. In the following discussion all data 
have been converted to particle sizes assuming a density of 2.4 g cm 3 . Whilst this may 
be a reasonable assumption for the density of individual grains, particles with high 
porosity may have a significantly lower bulk density. 

Firstly, to summarise the results here, there are two cumulative size 
distributions calculated from the features in the aerogel (see above and Table 6). The 
first is based on the Level 2 scan of surface features of the whole collector tray and 
the second is based on the higher resolution Level 3 scan of tracks (side view) in the 
aerogel blocks (15% of the total) removed during the PE period. The two thus use 
different data sets, the former is complete at large sizes (up to nearly mm impactor 
size) but loses sensitivity at intermediate sizes (10 pm impactor size), whilst the latter 
is never complete per se, but is sensitive to smaller impact features extending the 
impactor size range down to 0.5 pm. The Level 2 and 3 data sets are shown on the 
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same axes in Figure 23. A fit to the combined Level 2 and 3 results was made over the 
size range 1 to 200 (im in projectile size, and gave (dotted line on Figure 23): 
Cumulative size distribution per nr = 14100D P ' 0 ' 86±0 ' 01 . eqn (22) 

It can be seen on Figure 23 that at small sizes the aerogel data rolls off and is probably 
incomplete due to limited resolution. At large sizes a single fit no longer describes all 
the data (with the data falling below the fit result). 

These results can be compared to those previously reported by the Stardust PE 
team. In Horz et al., 2006, Fig. 4 shows the cumulative size distribution data. The fit 
that was given in that paper was based upon the craters measured in the aluminium 
foils carried by Stardust. Using the calibration in Kearsley et al. (2006), the crater 
widths (rim crest to rim crest) were converted to impactor diameter. This calibration 
used impacts of soda lime glass beads onto Stardust grade foils (similar in nature to 
those used here for the aerogel, with diameters in the size range 10 - 84 pm). Two 
sets of Stardust craters were used in obtaining the impactor size distribution. The first 
were called “large” craters (> 10 pm dia.) and represent the tray-wide sum of such 
craters on all foils, akin to the Level 2 track observations. The second category were 
small craters (<5 pm dia.) and were obtained from higher resolution SEM scans of 
selected sub-areas of the foils that were harvested and allocated during PE. After 
application of the calibration (i.e., assuming the impactors were spheres of density 2.4 
g cm" ) the small craters correspond to impactors of size 20 nm to 1pm (and thus are 
smaller than the projectiles used in the calibration) and the large craters to impactors 
of size 4 to 100 pm (roughly covering the calibration projectile size range). A single 
fit to the combined crater data over the impactor size range 0.05 to 100 pm yielded: 
Cumulative size distribution per m 2 = 125800 Dp‘ L72±0 ' 05 . eqn (23) 

The crater data are shown in Figure 23 (the fit is the solid line). Above a particle 
diameter of 10 pm the aerogel track and foil crater data overlap substantially. 
However, below this size there is a divergence that increases as particle size 
decreases. 

In Horz et al. (2006) the aerogel data that was used agrees at large particle 
sizes with that presented here, but differs below 10 pm. In the Level 3 scan data here 
all the extracted blocks have been used to obtain the cumulative size distribution. 
However, in Horz et al. (2006) only two blocks were used, namely blocks CO 12 and 
C023. As can be seen from Figure 22, these are the only two blocks with large 
numbers of tracks and the latter were dominated by small tracks and hence small 
particles. As was shown in Horz et al. (2006) (and reproduced here on Figure 23) the 
flux based on just these blocks is fully compatible with that from the foil crater data 
down to particle sizes of 1 or 2 pm. Below this, the aerogel data rolls off, indicating 
the limits of resolution have been reached. The reason why the track data from 
aerogel blocks CO 12 and C023 differs from that given by the Level 3 analysis based 
here on more blocks, is that the former appear to contain a cluster of impacts which 
has occurred on a spatial scale of cm across the collector and this is the only such 
cluster in the larger sample. Given that this cluster is dominated by small particles the 
cumulative size distribution is accordingly influenced at small sizes in the two 
analyses by the relative areas considered. As described in detail by Westphal et al. 
(2007) this clustering is poorly understood at present. 

The results obtained during the post flight analysis from the aerogel and foils 
can be compared to those obtained during the cometary encounter by the real time 
Dust Flux Monitor Instrument (DFMI - Tuzzolino et al. 2003). The data obtained by 
DFMI during the encounter were analysed in Tuzzolino et al. (2004) and in more 
detail in Green et al., 2004. Two points emerged. Firstly, the slope of the cumulative 
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size distribution varied as the spacecraft passed closer to the comet (within 3700 km, 
referred to as inner coma) and then travelled further away (outer coma). Secondly, 
above particle sizes of approximately 100 pm an excess of particles was reported 
which was not included in any fits. For comparison to the data here the DFMI data 
were converted from particle mass (m) to particle size (r) (with density 2400 kg m" 3 ) 
and are shown plotted on Figure 23 along with all the other data sets.. A fit to the 
DFMI data from 1 < r < 100 pm, yielded: 

Cumulative size distribution per nr = 1622000 D p ‘ 2 ' 55±0 ' 15 , eqn (24) 

(with regression coefficient -0.9785). As expected the power -2.55 is equal to the 
power of the equivalent mass distribution (-0.85) reported by Tuzzolino et al. (2004) 
and Green et al. (2004). In Figure 23 above 10 pm in particle size, DFMI, craters and 
aerogel tracks all yield similar results. The discrepancies between the data sets only 
emerge significantly at smaller particle sizes with extreme divergence between DFMI 
and the other data sets showing up at 1 pm. In the DFMI data the results for the 
smallest size particles were obtained from the most sensitive impact detector which 
had a surface area of 20 cm 2 . 

When comparing the size distribution measured during encounter and those 
measured after the return, the large difference in the power of the cumulative size 
distributions (Table 6) can be seen to lie in two effects. Firstly, the fits in the foil and 
aerogel analyses include data for impactors of size greater than order 100 pm. In the 
fits to the DFMI encounter data (Tuzzolino et al. 2004; Green et al. 2004) this part of 
the size distribution was excluded as it was considered to be due to an excess of large 
grains (similar to that reported by the Giotto encounter with comet Halley, e.g. 
McDonnell et al. 1991 and also partly as it was of low statistical significance, 7 
impacts). If the data at large sizes ( r > 100 pm) were included, the data point at the 
smallest sizes is then not compatible with the fit of a single power law. Accordingly if 
we remove the data point at the smallest size but include those at the largest sizes, we 
obtain: 

Cumulative size distribution per m 2 = 182000 D p ~* 74±015 , eqn (25) 

(with regression coefficient of -0.9787). There are thus two equally statistically valid 
fits (eqn.s 24 and 25) but with different slopes, dependent on which end of the size 
range is excluded. Based on comparison with the other data here, it is the DFMI 
datum at smallest sizes that appears anomalous, as even allowing for clustering events 
the other detectors do not reproduce such a high flux at small sizes and the DFMI data 
at this size are spread over a time period incompatible with a single cluster impact. In 
Green et al. (2004) the DFMI data were divided into time windows (corresponding to 
different stages of the encounter with the comet). The power of the size distribution 
was found to vary greatly depending on the time interval chosen, ranging from -0.9 to 
-2.25 near closest encounter, rising to -3.39 some 600+ seconds after the encounter 
when a burst of counts was detected. It should also be noted that the other active 
instrument during encounter CIDA (Cometary and Interstellar Dust Analyzer with 
very small active area) was designed to produce time of flight mass spectra from 
impacts of small particles (size > 0.1pm), and it observed a flux significantly smaller 
than that predicted by the other methods extrapolated to such scales (Kissel et al. 
2004); the single cumulative datum at approximately 1 um projectile size of CIDA 
(see Fig. 23) falls below most of the post-flight crater and track analysis by 
approximately a factor of 2 - 10. 

Taken together, the various DFMI, foil and aerogel data sets indicate that there 
was significant temporal variation in the particle flux and size distribution during 
encounter (DFMI) and that there were significant in-homogeneities in the spatial 
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distribution of the cumulative particle flux measured upon sample return (foil craters 
and aerogel tracks), particularly at smaller sizes, e.g. 1 pm scale and below. Obtaining 
a single flux or cumulative size distribution is thus difficult, as it will depend on the 
time period considered, the particle size range used and/or the location on the 
spacecraft of the detector. However, appropriate cumulative size distribution can be 
obtained if a set of constraints are applied. Namely: (1) Use the combined flux from 
the encounter and ignore any temporal variations. (2) Include the particle size range 
up to largest sizes (mm scale), i.e. do not treat large grains as an excess. The various 
data sets here that satisfy these constraints produce a cumulative size distribution that 
is compatible across all the individual data sets for sizes >10 pm and is best described 
by eqn. 13. (3) For sizes below this, discrete clusters of particles are observed and 
how many are included in the data set and normalised to what collecting area can 
greatly influence the reported cumulative size distribution. The range of slopes shown 
in Figure 23 at small sizes reflects variation in these parameters acting differently in 
the various data sets. The three fits shown in figure 23 illustrate this. A small number 
of clusters (potentially one) normalised to a large collection area yields a distribution 
with a shallow slope (eqn. 22). Data with a few clusters and a moderate area, yields an 
intermediate slope (which is well represented by the extrapolation of eqn. 13 to 
smaller sizes), whilst use of several clusters in a normalising area restricted to the 
cluster regions increases the slope to that of eqn. 23 or greater (e.g. eqn. 24). 
However, it should be noted that at the smallest sizes measured by DFMI the flux was 
not compatible with the other measurements, either in absolute number nor in the time 
distribution of the signals (which are incompatible with coincident hits inside a 
cluster). 

These results can be compared to the cumulative size distributions from 
previous cometary encounters. The power of the size distribution from 1 P/Halley was 
given as approximately -3.0 (McDonnell et al. 1991) or -(2. 6+0. 2) (Fulle et al. 2000), 
although this was found to depend on both particle size and distance from the nucleus. 
By contrast, the dust grain size distribution measured during encounter at 26P/Grigg- 
Skjellerup had a slope of -0.93 (McDonnell et al. 1993). Ground based observations 
of cometary dust mass distributions have also been made. For example, based on 
observations of comet Hale-Bopp combined with modelling and simulations, it is 
predicted that the small particle distribution in the size range 0. 1 to 20 pm had a size 
distribution slope of -3 (Levasseur-Regourd et al. 2007). The time-resolved Stardust 
data suggest considerable spatial and temporal variation in the size distribution of 
freshly liberated dust from a single comet, while the cumulative Wild 2 observations 
seem to suggest differences from comet to comet. 


5. Conclusions 

The Stardust mission successfully returned cometary dust grains captured in aerogel 
after being freshly emitted from comet P/Wild-2. Analysis of the tracks in the aerogel 
has permitted an estimate of the total cometary particle flux intercepted by the dust 
collector. The tracks in the aerogel divide into several categories depending on the 
variable composition and structure of the particles. There are solid grains which 
remain mostly intact during capture and more friable grains with either a less cohesive 
structure or a substantial volatile content. The latter may also contain smaller grains of 
well consolidated materials. 
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Laboratory impacts into aerogel have provided calibrations which were used 
to obtain Stardust cometary particle cumulative size distributions and provided 
samples as references for composition analysis teams. Study of 15% of the aerogel 
yielded 177 well identified tracks over 50 pm in length, which extrapolated to the 
whole cometary tray suggests that some 1180 tracks of this length or above are 
contained in the aerogel blocks. This is a rich harvest of cometary materials awaiting 
detailed analysis. 

The cumulative size distribution and flux obtained for Wild-2 shows non- 
uniform features in both short time and spatial regimes. For particle sizes greater than 
10 pm all the measurement methods used by Stardust produced similar results. 
However at smaller sizes, measurement of the particle size (or mass) distribution by 
different methods produced significantly different results that are not readily 
explained at this time; they may relate to the differing detection thresholds of the 
diverse methods or to spatial and temporal heterogeneities of the coma dust at the 
scale of individual detector surfaces e.g. stream effects in the coma and localised 
point sources (some of which may be close to the spacecraft, i.e., break up of larger 
particles after emission from the comet nucleus as discussed by Westphal et al. 2007). 
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Table 1: Properties of glass beads used in calibration work. Impact speeds are 
accurate to better than 1% 


Nominal Diameter 
±c (pm) 

Measured Diameter +c (pm) 
- see Kearsley et al., 2006 

Impact Speed 
(km s' 1 ) 

11.58 + 0.19 

9.8 + 2.2 

6.07 

35.0 + 0.8 

34.7+1.0 

5.99 

63.8 + 0.8 

64.1+2.8 

5.82 


Table 2: Average track dimensions measured in the aerogel. 


Nominal 
Diameter 
+c (pm) 

Number 

of 

tracks 

Track 

length 

(pm) 

Entrance 

hole 

diameter 

(pm) 

Maximum 

width 

(pm) 

Depth of 
maximum 
width 
(pm) 

Captured 

particle 

diameter 

(pm) 

Track volume 
(mm 3 ) 

11.58+0.19 

6 

2131+529 

84+35 

114+33 

238+130 

15+2 

0.010 + 0.007 

35.0 + 0.8 

5 

9309+1755 

178+22 

467+73 

2026+769 

34.8±1.8 

0.643 + 0.235 

63.8 + 0.8 

5 

17598+837 

379+76 

750+144 

1108+331 

66.0+2.4 

3.10+1.32 


Table 3: Mineral and basalt shots into Stardust grade aerogel in the laboratory. Impact 
speeds are accurate to better than 1% 

Projectile Source Impact Speed Grain size (pm ) 1 

(km s' 1 ) 


Basalt 

Natural History Museum 
(London). USGS sample 
NKT - 1G 

6.09 

1 - several hundred 
pm 

Pyrrhotite 

Natural History Museum 
(London) 

BM. 2005, M3 17 

5.85 

< 125 

Lizardite 

Natural History Museum 
(London). Confirmed by 
X-ray diffraction. 

5.94 

>53 

Olivine 

LLNL, San Carlos olivine 

5.85 

<45 


1 This was at launch, some minerals may break apart into a smaller size range during 
the shock of launch. 


Table 4: Cometary aerogel blocks used in this work. See Figure 8 for the 
layout of the blocks in the collector tray. 

Cometary aerogel trays used in this work 

C009, C012, C013, C023, C027, C029, C038, C048, C049, C052, C054, C086, C091, 
C092, C101, C102, 018, 026, 027, 028 


Table 5. Track categories vs. track length (where % do not sum to 100 this is due to 


round- off error). 


Track length range 

% Type A 

% Type B 

% Type C 

All 

65 

33 

2 

<100 pm 

100 

0 

0 

100 - 500 pm 

70 

27 

3 

500 - 1000 pm 

49 

49 

3 

> 1000 pm 

51 

49 

0 
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Table 6: Fits to cumulative size distribution, of the form y = ax’, with x in pm. 


Data set 

Impactor 
range (pm) 

size a 

b 

Aerogel Level 2, eqn. 13 

10 - 2000 

47100 

-1.22+0.06 

Aerogel Level 3, eqn 21. 

1 - 100 

11900 

-0.76+0.01 

Aerogel (all), eqn. 22 

1-200 

14100 

-0.86+0.01 

Foil Craters, eqn. 23 here 
and Horz et al., 2006. 

0.05 - 100 

125800 

-1.72+0.05 

DFMI, Tuzzolino et al., 
2004 and eqn. 24. 

1-200 

1622000 

-2.55+0.15 

DFMI (excluding small 
sizes and including large 
szies) eqn. 25 

5-600 

182000 

-1.74+0.05 
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Figures 


Figure 1 . Schematic of typical “carrot” shaped track observed in laboratory impacts of 
glass beads in aerogel. Key dimensions are labelled. (EHD = entrance hole diameter). 
Note that the width at track end is not always the same as the mean diameter of the 
captured particle. Impact direction was from left. 

Figure 2. For the calibration shots in aerogel: (a) Average track length vs. original 
particle diameter, (b) Average entrance hole diameter vs. original particle diameter, 
(c) Average greatest track width vs. original particle diameter. Parameterisations of 
the fit curves are given on the figures (and explained in the main text as eqn.’s 1, 3, 5 
and 7) along with r, the regression coefficient of the fit. 

Figure 3. (a) Top and side views of a frustum shape, (b) Side view of two frustums, 
which approximate the carrot track shape shown in Figure 1 . 

Figure 4. The estimated average track volume vs. impact kinetic energy is shown for 
the calibration shots in aerogel. The parameterisation of the fit curve (eqn. 11, see 
main text) is given on the figure along with r, the regression coefficient of the fit. 

Figure 5. Tracks made by impact of: (a) pyrrhotite grains in aerogel at 5.85 km s' 1 ; the 
largest track just extends off the bottom of the image so the captured grain is not 
visible, (b) lizardite grain in aerogel at 5.94 km s' 1 ; the incident grain has broken 
up/disaggregated into many fine fragments which are too small to be visible at this 
scale but line the walls of the cavity in the aerogel. In both cases the impact direction 
was from the top. 

Figure 6. Top down view of the entrance hole in aerogel made by the impact of a soda 
lime glass bead (63.8 pm diameter) at 5.82 km s' 1 in the laboratory. 

Figure 7. Layout of aerogel blocks in cometary aerogel collector tray. Each block is 
labelled with its Stardust designator Cnnn, where “nnn” is an integer in the range 1 - 
132. 

Figure 8. Example entrance holes in Stardust cometary aerogel (track C027-T4). The 
edge of the hole in the top surface is shown indicated by an arrow (A), the edge of the 
region of widest extent in the subsurface is also shown indicated by an arrow (B). 

Figure 9. Size frequency of features seen on Level 2 scan. In addition, there are 4 
features with size greater than 4 mm. 

Figure 10. Distribution of Level 2 features across the aerogel collector. Feature 
diameter is represented by symbol size: ■ is 100 - 500 pm, • is 500 - 1000 pm, • is 

1000 - 5000 pm and * is 5000 - 10000 pm. 

Figure 11. Cumulative particle size distribution for the cometary aerogel collector; 
based on an optical (Level 2) scan of the whole aerogel surface whilst it was still in 
the collector tray (before any individual aerogel cells were removed from analysis). 
The fit curve shown is described in the main text (eqn. 13). 
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Figure 12. Examples of Type A impacts from Stardust cometary aerogel tray (Level 3 
images), (a) Typical track (C038-T7). (b) Captured terminal particle (bottom) at end 
of track (from C086-T2). In both cases, impact direction was from the top. 

Figure 13. Example of Type B impact (C054-T3) from Stardust cometary aerogel 
tray (Level 3 images). Impact direction was from the top. 

Figure 14. Example of Type C track (C052-T6) from Stardust cometary aerogel tray 
(Level 3 images). Impact direction was from the top. 

Figure 15. The % of tracks that are Type A, defined as 100 x EA / E(A+B) vs. track 
length (bin width 100 pm). For the longest tracks the quantization of 0, 66 or 100% at 
long track lengths is due to small numbers of tracks per bin. 

Figure 16. Frequency of the number of terminal tracks (“styli”) for (a) Type A tracks 
and (b) Type B tracks. 

Figure 17. Stardust cometary aerogel track widths vs. lengths (Level 3 scan). Track 
types are as indicated. See main text for details of the boundary lines superimposed on 
the data. 

Figure 18. Ratio of track width/length vs. track length. Track types are as shown. See 
main text for details of the boundary lines superimposed on the data. 

Figure 19. Track volume (V) vs. track length (L). for tracks measured in the Stardust 
cometary areogel: (a) Type A tracks, (b) Type B, (c) all data. The fit curves shown in 
(a) and (b) are eqn.s 16, 17, 18 as described in the main text. The boundary lines 
shown in (c) are eqn.s 19 and 20 as described in the main text. 

Figure 20. Estimated particle size vs. track length for tracks measured in the Stardust 
cometary aerogel. Particle size was obtained from eqn. 12 in the main text. 

Figure 21. Cumulative particle size distribution based on a high resolution optical 
(Level 3) scan of the 20 aerogel blocks extracted from the Stardust cometary dust 
collector tray. The fit curve shown is described (eqn. 21) in the main text. 

Figure 22. Spatial distribution of tracks in Level 3 scan of extracted aerogel cells, 
aerogel. Extracted cells are shown shaded, the number in those cells is the observed 
number of tracks. Note the large numbers of tracks in cells 12 and 23 are 
predominantly very short tracks. 

Figure 23. Cumulative particle size distribution based from all data sets: aerogel level 
2 and 3 scans (data presented herein), aerogel cells 12 and 23 (Horz et al., 2006), 
craters in foils (Horz et al., 2006) and the Dust Flux Monitor Instrument (DFMI, 
Tuzzolino et al., 2004). The parameterisations of the fit curves shown are given in 
Table 6. 
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Figure 1 . Schematic of typical “carrot” shaped track observed in laboratory impacts 
of glass beads in aerogel. Key dimensions are labelled. (EHD = entrance hole 
diameter). Note that the width at track end is not always the same as the mean 
diameter of the captured particle. Impact direction was from left. 
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Figure 2. For the calibration shots in aerogel: (a) Average track length vs. original 
particle diameter, (b) Average entrance hole diameter vs. original particle diameter. 
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(c) Average greatest track width vs. original particle diameter. Parameterisations of 
the fit curves are given on the figures (and explained in the main text as eqn.’s 1, 3, 5 
and 7) along with r, the regression coefficient of the fit. 
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(b) 



Figure 3. (a) Top and side views of a frustum shape, (b) Side view of two frustums, 
which approximate the carrot track shape shown in Figure 1. 
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Impact Kinetic Energy (J) 


Figure 4. The estimated average track volume vs. impact kinetic energy is shown for 
the calibration shots in aerogel. The parameterisation of the fit curve (eqn. 11, see 
main text) is given on the figure along with r, the regression coefficient of the fit. 
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Figure 5. Tracks made by impact of: (a) pyrrhotite grains in aerogel at 5.85 km s' 1 ; the 
largest track just extends off the bottom of the image so the captured grain is not 
visible, (b) lizardite grain in aerogel at 5.94 km s' 1 ; the incident grain has broken 
up/disaggregated into many fine fragments which are too small to be visible at this 
scale but line the walls of the cavity in the aerogel. In both cases the impact direction 
was from the top. 
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Figure 6. Top down view of the entrance hole in aerogel made by the impact of a 
soda lime glass bead (63.8 pm diameter) at 5.82 km s' 1 in the laboratory. 
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Figure 7. Layout of aerogel blocks in cometary aerogel collector tray. Each block is 
labelled with its Stardust designator Cnnn, where “nnn” is an integer in the range 1 - 
132. Shaded blocks were used in this study. 


32 


















































Figure 8. Example entrance holes in Stardust cometary aerogel (track C027-T4). The 
edge of the hole in the top surface is shown arrowed (A), the edge of the region of 
widest extent in the subsurface is also shown arrowed (B). 
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Figure 9. Size frequency of features seen on Level 2 scan. In addition, there are 4 
features with size greater than 4 mm. 
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Figure 10. Distribution of Level 2 features across the aerogel collector. Feature 
diameter is represented by symbol size: ■ is 100 - 500 pm, • is 500 - 1000 pm, • is 

1000 - 5000 pm and * is 5000 - 10000 pm. 
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Figure 11. Cumulative particle size distribution for the cometary aerogel collector; 
based on an optical (Level 2) scan of the whole aerogel surface whilst it was still in 
the collector tray (before any individual aerogel cells were removed from analysis). 
The fit curve shown is described in the main text (eqn. 13). 
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Figure 12. Examples of Type A impacts from Stardust cometary aerogel tray (Level 3 
images), (a) Typical track (C038-T7). (b) Captured terminal particle (bottom) at end 
of track (from C086-T2). In both cases, impact direction was from the top. 
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Figure 13. Example of Type B impact (C054_T3) from Stardust cometary aerogel 
tray (Level 3 images). Impact direction was from the top. 
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Figure 14. Example of Type C track (C052_T6) from Stardust cometary aerogel tray 
(Level 3 images). Impact direction was from the top. 
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Figure 15. The % of tracks that are Type A, defined as 100 x 2 A / E(A+B) vs. track 
length (bin width 100 pm). For the longest tracks the quantization of 0, 66 or 100% at 
long track lengths is due to small numbers of tracks per bin. 
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Figure 16. Frequency of the number of terminal tracks (“styli”) for (a) Type A tracks 
and (b) Type B tracks. 
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Figure 17. Stardust cometary aerogel track widths vs. lengths (Level 3 scan). Track 
types are as indicated. See main text for details of the boundary lines superimposed on 
the data. 
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Figure 18. Ratio of track width/length vs. track length. Track types are as shown. See 
main text for details of the boundary lines superimposed on the data. 
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(a) and (b) are eqn.s 16, 17, 18 as described in the main text. The boundary lines 
shown in (c) are eqn.s 19 and 20 as described in the main text. 
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Figure 20. Estimated particle size vs. track length for tracks measured in the Stardust 
cometary aerogel. Particle size was obtained from eqn. 12 in the main text. 
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Figure 21. Cumulative particle size distribution based on a high resolution optical 
(Level 3) scan of the 20 aerogel blocks extracted from the Stardust cometary dust 
collector tray. The fit curve shown is described (eqn. 21) in the main text. 
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Figure 22. Spatial distribution of tracks in Level 3 scan of extracted aerogel cells, 
aerogel. Extracted cells are shown shaded, the number in those cells is the observed 
number of tracks. Note the large numbers of tracks in cells 12 and 23 are 
predominantly very short tracks. 
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Figure 23. Cumulative particle size distribution based from all data sets: aerogel level 
2 and 3 scans (data presented herein), aerogel cells 12 and 23 (Horz et al., 2006), 
craters in foils (Horz et al., 2006) and the Dust Flux Monitor Instrument (DFMI, 
Tuzzolino et al., 2004). The parameterisations of the fit curves shown are given in 
Table 6. 
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